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IDEAL DECOMPOSITIONS AND COMPUTATION OF TENSOR 

NORMAL FORMS 

BERND FIEDLER 


Abstract. Symmetry properties of r-times covariant tensors T can be de¬ 
scribed by certain linear subspaces W of the group ring K[5j.] of a symmet¬ 
ric group Sr- If for a class of tensors T such a, W is known, the elements of 
the orthogonal subspace of W within the dual space IK[5r.]* of K[5r] yield 
linear identities needed for a treatment of the term combination problem for 
the coordinates of the T. We give the structure of these W for every situa¬ 
tion which appears in symbolic tensor calculations by computer. Characterizing 
idempotents of such W can be determined by means of an ideal decomposition 
algorithm which works in every semisimple ring up to an isomorphism. Fur¬ 
thermore, we use tools such as the Littlewood-Richardson rule, plethysms and 
discrete Fourier transforms for Sr to increase the efficience of calculations. All 
described methods were implemented in a Mathematica package called PERMS. 


1. The Term Combination Problem eor Tensors 


The use of computer algebra systems for symbolic calculations with tensor ex¬ 
pressions is very important in differential geometry, tensor analysis and general 
relativity theory. The investigations of this paper]] are motivated by the follow¬ 
ing term combination problem or normal form problem which occurs within such 
calculations. 

Let us consider real or complex linear combinations 

n 

(1.1) r = , Q!ieM,C 

i=l 

of expressions T(p which are formed from the coordinates of certain tensors A, R, 
C, ... by multiplication and, possibly, contractions of some pairs of indices. An 
example of such an expression is 


( 1 - 2 ) . 

Iii(Q) we use Einstein’s summation convention. Further we assume that each of 
the numbers of A, B,C,... is constant if we run through the set of the Up). Now 
we aim to carry out symbolic calculations with expressions of the type (p..l|), (|1.2|) 
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according to the rules of the Ricci calculus. We assume that there is a metric 
tensor g which allows us to raise or lower indices, for instance 


rpjkb _ bcrpjk T' ^ _ n T 

■^2 y i c "! ^ ia k yac^i k * 


If now the tensors A, B, C. have symmetries and/or fulhl linear identities^, 
then there exist relations between the T(j) in (|1 . 1|) . (We restrict us to linear rela¬ 
tions.) Thus the problem arises to detect such relations in sums ( |1.1|) , generated 
by symbolic calculations, and to reduce (|1.1|) to linear combinations of linearly 
independent r(j) {normal forms). 

It is well-known that the representation theory of symmetric groups Sr yields 
powerful tools to treat this term combination problem. The connection between 
tensors and the representation theory of Sr has been considered already in books 
of J.A. Schouten (1924), H. Weyl (1939) and H. Boerner 0 (1955). In the 
1940s Littlewood has developed and used tools such as the Littlewood-Richardson 
rule and plethysms for the investigation of tensors (see S.A. Fulling et ah |T4 
(references) and D.E. Littlewood (appendix)). 

Applying the same methods. Fulling, King, Wybourne and Cummins [M have 
calculated large lists of normal form terms of polynomials of the Riemann curvature 
tensor and its covariant derivatives (by means of the program package Schur [P^). 
In their paper they formulated the following steps to solve the above term 
combination problem for tensors: 


(a) Generate the space spanned by the set of homogeneous monomials of a 
dehnite ’order’ or ’degree of homogeneity’ formed from the coordinates of 
tensors of relevance by multiplication and index-pair contraction. 

(b) Construct a basis of this space {normal forms). 

(c) Present an algorithm for expressing an arbitrary element of the space in 
terms of the basis. 


Our present paper yields a method to solve (b) and (c) for arbitrary tensors. 


2. Tensors and the Group Ring of a Symmetric Group 

We make use of the following connection between r-times covariant tensors 
T E TrV over a hnite-dimensional K-vector space V and elements of the group 
ring K[iSr] of a symmetric group Sr over a held K = M or IK = C. 


Definition 2.1. Any tensor T G 7/R and any r-tuple b := {vi,... ,Vr) G W of r 
vectors from V induce a function Tb : 5,. — > IK according to the rule 

( 2 . 1 ) Ti){p) . T{Vp{\p ■ ■ ■ 1 '^p(r)) ) P £ Sr . 

We identify this function with the group ring element := '^^{p) P £ IK[iSr]. 

^For instance, the Riemann tensor has the symmetry Rijki = —Rjiu = —Rijik = RkUj and 
fulfils Rijki + Rikij + Riijk = 0. See P. Gunther [|^ for a large collection of identities for the 
curvature tensor, Weyl tensor, etc. and covariant derivatives of these tensors. 
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We allow the linear dependence of the Vi and repetitions of vectors in the above 
r-tnple b. Obviously, two tensors S,T E%V fulhl S' = T iff S';, = Tf, for all h . 

We try to describe symmetry properties of tensors with the help of the Ti, by 
the following principle: If a ’’class” of tensors with a certain symmetry property is 
given, then we search such a linear subspace W C ]K[5r] that contains all of the 
tensors from the ’’class” being considered. The linear identities which characterize 
W yield then identities for the coordinates of the T which can be used in the 
treatment of the term combination problem. 

Every a = J2p£Sr P ^ acts as so-called symmetry operator a :T ^ aT 

on tensors T G if we dehne 

(2.2) {aT){vi,...,Vr) := ^ a(p) T(np(i),..., Upp)) , Vi^V. 

p&Sr 

We denote by ’*’ the mapping *: a = J2p&Sr^^P^P 

Furthermore, if p E Sr and b = {vi,... ,Vr) G E”, then we denote by pb the r-tuple 
pb := (np(i),..., Vp(r))- Many of our calculations are based on 


Lemma 2.2. ^ If a = Q G T G TrV , p,q G Sr and b = 

{vi,... ,Vr) G W, then we have 


(2.3) 

(2.4) 

(2.5) 

( 2 . 6 ) 


Tb{p o q) 

q{pb) 

(aT), 

n 


{qT)b{p) = Tpb{q) 

{p o q)b 

Tb-a* 

P-Tpb. 


The following symmetry concepts are used for tensors. (See B. Fiedler 0. 
Sec.III.2] and B. Fiedler |^. See also R. Merris |^, pp.151,153,157], H. Boerner 
10, p.l27], G. Eisenreich 0, p.601].) 


Definition 2.3. (a) Let r C ]K[iSr] be a right ideal of ]K[5r] for which an a G r 

and a T G TrV exist such that aT ^ 0. Then the tensor set T( := {aT \ a G 
r, T G TrV} is called the symmetry class of tensors dehned by r. 

(b) Let ai,..., a„ G ]K[iSr] be a hnite set of group ring elements. We say that 
a tensor T G TrV possesses a symmetry dehned by Oi,..., if T satishes 
the linear equation system aiT = 0, (z = 1,..., n). 

If e is a generating idempotent of a right ideal r = e ■ ]K[5r] that dehnes a 
symmetry class T then % fulhls T^ = {eT \ T E TrV} and a tensor T E TV 
belongs to T( iff eT = T. (See H. Boerner 0, p.l27] or B. Fiedler []^ Sec.III.2.1].) 

Now it can be shown that all T, of tensors T which have a symmetry (a) or (b) 
lie in a certain left ideal of ]K[i5r]. 


^See B. Fiedler |10, Sec.III.1] and B. Fiedler [ p^ . 
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Proposition 2.4. ^ Let e G K[i5r] be an idempotent. Then a T E %V fulfils the 
condition eT = T iff G i := K[5r] • e* for all h G , i.e. all Tf, of T lie in the 
left ideal [ generated by e*. 

Proposition 2.5. [| Let ai,. .. ,am G K[5,.] be given group ring elements. AT E 
%V satisfies a system of linear identities OiT = 0 , {i = 1,... ,m), iff 
Tb E { ■.= {x E K[5r] \ X ■ a* = 0 , i = 1,..., m} for all b E , i.e. all T^ ofT lie 
in the left annihilator ideal [ of the set {af ..., a'^}. 


The proofs follow easily from (|2.5|) . A further result is 


Proposition 2.6. Q //dimP > r, then every left ideal 1 C ]K[iSr] fulfils 1 = 
£K{Tb \ T E Ti* , b E V^}. (Here Ck denotes the forming of the linear closure.) 

If dim V < r, then the Tf, of the tensors from Ti* will span only a linear subspace 
of I in general. 

In the case of tensors T with index contractions the role of the Tf, is played by 
certain sums XlbeSt lb Tb, which we now define. 

Definition 2.7. Let g E T 2 V be a fundamental tensor with arbitrary signature on 
V and B = {ni,..., n^} be an orthonormal basis of V with respect to g. Further 
let r, I be integers with 2 < 21 < r and bo = {v 2 i+i,... ,Vr) E be a fixed 

(r — 2/)-tuple of vectors from B. Then we denote by lBf,g the set [] of r-tuples of 
basis vectors 

:= {{Wi,Wi,W2,W2, • • . ,tCf,tCf,n2Z+i, ...,Vr)EB'^\ {wi,... ,wi) E B’^} . 

Moreover, we set jb ■= n!=i E {1, —1} for every b E ^bo- 


Proposition 2.8. ^ Let T E %V be a tensor and g E T 2 V be a fundamental 
tensor. We determine all tensor coordinates with respect to an orthonormal basis 
B = {rii,..., Ud} of V. Let bo = + 1 ,..., Ui^) E B^~‘^’‘ be a fixed (r — 2l)-tuple of 

basis vectors. Then 


(2.7) 


peSr 


n 32 
31 32 .. 


■■ 31 

jl *2i + l...*r 


P = 


be^bo 


Due to Prop. || the 'ffb&<Bb 1^ objects which contain information about 

the tensor coordinates of T with I index-pair contractions. In the case of tensors T 


12 


'^See B. Fiedler 
^See B. Fiedler 

U 


®See B. Fiedler 
'^We set bo := 
r = 2l>0. 


or B. Fiedler |K, Prop. III.2.5, III.3.1, III.3.4]. 
or B. Fiedler jw, Prop. III.3.3, III.3.4]. 
or B. Fiedler |1C, Prop. III.2.6]. 

and $0 := { {wi,Wi,W 2 , W 2 ,..., w/, wj) E | (wi ,... ,wi) E B’^} in the case 


See B. Fiedler |12 or B. Fiedler Prop. III.3.31]. 
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with index contractions we search for linear subspaces W C ]K[i5r] which contain 

all ^ ^0- 

The left ideals 1 from the Propositions I, I and 1^ are the simplest examples of 
linear subspaces W describing tensor symmetries. Before we give such subspaces W 
in more complicated cases, we will explain how they can be used in the treatment 
of the term combination problem for tensors.|^ 


3. The Treatment of the Term Combination Problem 


The term combination problem from Sec. |I] can be reformulated as the problem 
to hnd all linear identities between the summands of given tensor expressions 


(3.1) 

(3.2) 




T, 


l2i + l...*r 


^ ^ ^p(l) •••^p(r) 

p&V 

pGV 


Cp G IK, "P C or 


(PT) 


n J 2 

31 32 ■■ 


■■ 31 

jl i2l + l-"ir 


G IK , P C Sr 


where T G p* is a tensor from a given symmetry class dehned by a left ideal 1 
(or a right ideal r = [*). We assume that ( p.lf ) and are results of symbolic 
computer calculations. P is a subset of permutations, which is determined by the 
concrete form of the given expressions (|3.1|) or (|3.2|) . 

Let VP C 1 be a linear subspace which contains all p or Tfc P of T. If we 

consider the orthogonal subspace := {x G K[iSr]* | Vtc G VP : (t, u;) = 0} of 
VP, then every x G VP^ yields a linear identity for the coordinates of T since 

(3.3) 0 = {x,Tb) = '^XpTbip) = 

p€Sr p£Sr 


or 


0 = 


{x, IbTb) = Y^bTb{p)xp = 


be^bo 


bs^bf) 

p€Sr 


p&Sr 


X, (PT)P\. 


■■ 31 

jl i2l + l-"ir 


where Xp := {x,p), p G Sr- (The last steps are correct if all b occurring in (|3.3|) 
are r-tuples of basis vectors of V.) Every identity (|3.3| ) can be used to eliminate 
certain summands in ( |3.1|) , ( p.2|) . If VP is spanned by all p, or J2bG<Bb IbTb of 

the tensors considered, then VP^ contains all linear identities which are possible 
between summands of expressions (|3T|) , (|3.2|) (compare Prop. P). 


If a basis {hi,..., /i^} of VP is known, then the coefficients Xp of the x G VP^ 
can be obtained from the linear equation system 


(3.4) (t, hi) = = 0 = 

pSSr 

Thus an important goal is to hnd such a basis (hi,..., h^} of VP. (An efficient 
algorithm for that purpose is given in Prop. 0 .) 

Note that ( |3.4| ) is a very large system with a (h x r!)-coefficient matrix, k = 
dim VP. However, since we only need identities to reduce sums m , (|3^), we can 


®See B. Fiedler |12 or B. Fiedler ||T^, Sec. III.l, III.4.1]. 
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restrict us to solutions of (|3.4|) which fulhl Xp = 0 for p G \ V. This reduces 
the number of columns to \V\. Furthermore, every of our spaces VF is a linear 
subspace of a left ideal 1 = ]K[iSr] • e. A decomposition e = Ci + ... + of the 
generating idempotent into pairwise orthogonal, primitive idempotents e* induces 
a decomposition W = Wi © ... © Wm with Wi C ]K[5r] • e* and a decomposition of 
the tensors T E T\*\ T = e]T + ... + e*^T. Then we can transform (|0|), (|3^ ) into 
expressions formed from the e*T, for instance 


(3.5) (pD 


r, 


*2i + l- 


EE' 

i=l pGR 


(P«T)) 


31 

31 32 ■■ 


■■ 3l 

3l *2i+l---V ’ 


and use the smaller equation systems of the smaller spaces Wi to determine linear 
identities for the coordinates of the e*T. Note, however, that a reduction of expres¬ 
sions such as (|3.5|) by means of identities of the Wi leads to a linear combination 
of coordinates of the e*T which cannot be reckoned back into a linear combination 
of the coordinates of T in general. 


4. The Algorithm for Ideal Decompositions 

Two problems became visible up to now: We need methods 

(a) to determine generating idempotents e for left/right ideals of ]K[iSr] for 
which such idempotents are unknown (such as in Prop. 

(b) to decompose a given idempotent e G ]K[5r] into a sum e = ei + ... + of 
pairwise orthogonal primitive idempotents e*. 

We developed an algorithm which solves these problems. First versions of this 
algorithm were designed for ]K[5r] (see B. Fiedler [|, ||]). But it turned out that 
this algorithm works even in an arbitrary semisimple ring 91 which fulhls: 

(A) We know explicitly a decomposition 

m m 

(4.1) 9l = 09I-i/, or 91 = 01/,-91 

i=l i=l 

of the full ring 91 into minimal left or right ideals generated by known 
primitive idempotents i/,. Pairwise orthogonality of the y, is not required. 

(B) We know in 94 a method to construct explicitly a solution x G 94 for every 
equation 

(4.2) e ■ a ■ X ■ e = e or e ■ x ■ a ■ e = e , 

where e G 94 is a primitive idempotent and a G 94 is a ring element with 
e-OT^Oora-eT^O, respectively. 

We describe now a version (L) of the algorithm for left ideals of a semisimple ring 
94. Obviously an analogous version (R) for right ideals can be formulated, too. 
(See B. Fiedler |^, Sec.1.2] and B. Fiedler p, ^.) 

A frequent step of the algorithm is the construction of a generating idempotent 
for a left ideal 1 = 94 ■ e • a, where e G 94 is a primitive idempotent and a G 94 is a 
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ring element with e ■ a ^ 0. This is possible by 


Proposition 4.1. For the above e, a G 91 there exists such an x G 91 that 
(4.3) e ■ a ■ X ■ e = e . 

Moreover, the ring element e' := x ■ e - a formed with this x is an idempotent which 
generates the minimal left ideal 91 • e • a. 

A second construction orthogonalizes idempotents. Let I = 91 • e and 1 = 91 • e 
be two left ideals with known generating idempotents e and e. Assume that e is 
primitive and e • e 7 ^ e, i.e. I ^ 1. Then the sum 1 + 1 = 1 © 1 is direct since the 
minimality of I yields 1 fl 1 = {0}. Now we search for new generating idempotents 
/ and / of [ and I which fulhl f ■ f = / ■ / = 0 . 


Theorem 4.2. Q The above orthogonalization problem can be solved in 2 steps: 

(i) We can find such a ring element x G 91 that 

(4.4) e • (1 — e) • X • e = e . 

If we use this x to form f := (1 — e) ■ x ■ e, then f is a generating idempotent 
of 1 which satisfies e ■ / = 0. 

(ii) For a result f of step (i) there exists an x G 91 such that 

(4.5) /•(l-e)-x-/ = /. 

If we make use of x to form f := e — {1 — e) ■ x ■ f ■ e, then f is a generating 
idempotent of I which fulfils /■/ = /•/ = 0. 


Now we can describe our Algorithm (L), which allows us to decompose every 
left ideal I = 91-a with known generating element a 7 ^ 0 into a direct sum 1 = 0^^ h 
of minimal left ideals h explicitly (see B. Fiedler [|T^, Sec.1.2] and B. Fiedler P, |^). 
A multiplication of ( |4.1|) by a yields a sum 

m 

(4.6) 1 = 91-a= ^91-2/i-a 

i=l 

Vi-a^O 


of minimal left ideals for 1 which however is not direct in general. Now we can 
carry out the following steps: 


(1) The first summand in ( [4.6|) is a minimal left ideal. We denote it by b and 
we determine a generating idempotent ei of b by means of Prop. 0 . 


^^See B. Fiedler Prop. 1.2.1]. Compare B. Fiedler ||, ||]. 

^^B. Fiedler Thm. 1.2.4] and B. Fiedler ]^, |j. The proof uses the fact that the set 
{e — x-e + e- x- e\ x£ 91} is the set of all generating idempotents of 1 = 91 • e (compare D.S. 
Passman |^, p.l37]). 
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(2) We search for the first minimal left ideal IK ■ i/j ■ a in (|4.6|) which is not 
contained in h that means for which 


(4.7) 


yi - a-ei 7 ^ pi 


We denote it by h and we construct a generating idempotent 62 of I 2 by 
means of Prop. Since [2 is minimal and h ^ b , we obtain li fl I 2 = {0}. 
Thus b and b form a direct sum b := b © b- Because 62 is primitive and 
62-61 7 ^ 62 , we can determine new, orthogonal, generating idempotents 
fl, /2 of b 1 b by means of Theorem 0 . Then /2 := /i + /2 is a generating 
idempotent of b- 

(3) Now we search for the next minimal left ideal ■ a in ([4-61) which is not 
contained in b that means for which 


(4.8) 


Vi ' Oi ' f2 7 ^ Vi 


We denote it by b- We construct a primitive generating idempotent 63 of [3 
and pass over to new orthogonal idempotents / 2 , fs instead of / 2 , 63 . This 
leads us to the left ideal b := b © Is which has the generating idempotent 
/s := f2 + /s- 

(4) We continue this procedure until we have processed all left ideals in ( |4.6|) . 
The result is a left ideal b := b © • • • © b and a generating idempotent /„ 
of b- 

Obviously, b © 1 since every left ideal U is a summand in ( |4.6|) . Furthermore, every 
summand 94 ■ ■ a from ( |4.6| ) which had been considered before we had reached 

the ideal b is contained in b-i = b © • • • © b-i © b- All other summands 94 ■ ■ a 

of (|4.6|) lie in b. This leads to 1 C b and 1 = b- 

According to Theorem 0 (ii), every generating idempotent fk of b can be written 
sis fk = {1 — Zk) ■ fk with a Zfc G 94 which we have already determined to carry out 
the orthogonalization (cfc+i, fk) 1 —(A+i, fk)- Thus we can write 


In 


fn—1 © fn — (1 Zn—l) ' fn—1 © fn — (1 l) ' {In—2 © fn—l) © fn 

(1 Zn—l) ' (1 Zn— 2 ) ' In—2 © (1 Zn—l) ' fn—1 © fn 


n-1 

( 4 . 9 ) = '^{l - Zn-l) ■ {1 - Zn- 2 ) ■ ■■■ ■ {'^ - Zk) ■ fk + fn- 

k=l 

Formula (|4.9|) presents a decomposition fn = Ylk=i ^k of fn into summands which 
fulfil /ifc := (1 - Zn-i) ■ (1 - Zn-2) {I - Zk) ■ fk e b and hn := fn © b- 

Thus, fn = YIk=i decomposition of fn which corresponds to 1 = 0^=1 b 

and the hk are pairwise orthogonal generating idempotents of the b- 

Obviously, the algorithm (L) solves the above problem (b). Furthermore, the 
algorithm (L) can be extended to left ideals which are non-direct sums 1 = 
tti by applying its steps to the summands of 1 = Yli=i XlJLi ^' n*. Likewise, 
we can construct generating idempotents and decompositions for right ideals r = 
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^^^=1 Oi ■ 5^ by the algorithm version (R). If our left/right ideals are intersections 
1 = nf=i 9^ • e, or r = nf=i e* ■ of left/right ideals (e* idempotents), then their 
right/left annihilator ideals are Arii) = or ^^(r) = 

respectively. In this case we can construct a generating idempotent e of Ar{\) or 
Ai{t) by (R) or (L), respectively, and form e' := 1 — e to obtain a generating 
idempotent e' of I or r. Thus our algorithms solve problem (a) for non-direct sums 
or intersections of left/right ideals. See B.Fiedler Sec.1.4] or B. Fiedler ^ for 

further details. 

5. Completions of the Decomposition Algorithms 

The basic assumptions (A) and (B). Actual decomposition constructions 
can only be carried out by the algorithms (L) or (R) if our semisimple ring fH fulhls 
the above assumtions (A) and (B). This is the case for 

(1) the group ring = ]K[5r] of a symmetric group, 

(2) a ring IR = that is an outer direct product of full (n, x rij)- 

matrix rings over skew helds §*, 

(3) all semisimple rings IR for which an isomorphism D : IR —IR' = (S)™ i S/*^*^* 
onto a ring IR' of the second type is explicitly known. 

(See B. Fiedler |^, Sec.1.3].) According to Wedderburn’s Theorem every semisim¬ 
ple ring is isomorphic to a ring of type 2. Thus statement 2 means that the 
decomposition algorithms (L) and (R) work in every semisimple ring up to an 
isomorphism. 

Let us consider fR = We denote by Cki G a matrix in 

which exactly the element located in the fc-th row and the j-th column is equal 
to 1 G whereas all other elements vanish. Then a decomposition of 93 into 
minimal left/right ideals is given by the decompositions S”*^*^* = 0 JLi ■ Cjj 

and = 0"^^ (7^.^. . Qf matrix rings into minimal left/right ideals. 

Furthermore there exists a very fast procedure to solve ( [4.2|) in 93. 

Since a primitive idempotent e G 93 has only one non-vanishing block matrix 
(i.e. e = (0,..., A,..., 0), where E G is also a primitive idempotent), an 

equation such as e - a - x ■ e = e leads to a single matrix equation E ■ A - X ■ E = E. 
Moreover, E can be written as E = ■ h with row vectors f,hE , where /i is a 

non-vanishing row of E. If we set m = h - A and determine non-vanishing elements 
m^o, fko of /; then X := is a solution of E ■ A ■ X ■ E = E, 

which yields a solution x = (0,..., A,..., 0) of e ■ a • x ■ e = e (B. Fiedler |Ty, Sec. 
1.1.2, 1.3.2]). Obviously, this procedure will run very fast on a computer. 

For ]K[5r] the well-known decomposition of ]K[iSr.] into minimal left/right ideals 
by means of Young symmetrizers guarantees (A). See B. Fiedler p, P| for (B). 

If for a semisimple ring 93 an above isomorphism D : 93 —93' is known (and 
practicable on a computer), then every ideal decomposition problem for 93 can be 
transferred to 93' and treated there by the algorithms (L) and (R). 

If 93 = C[G] is the group ring of a hnite group G, then we have Sj = C for all i and 
the isomorphism D is called a discrete Fourier transform for G. Explicit algorithms 
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for such Fourier transforms are known at least for abelian groups, solvable groups, 
supersolvable groups and symmetric groups (see M. Clausen und U. Baum 0). 


Discrete Fourier transforms. In group rings = C[G] of large finite groups G 
(such as G = Sr, r > 8), even a single product a ■ b, a,b E C[G] can lead to high 
costs in time and computer memory (see B. Fiedler [^, Sec. 1.1.3, 1.5.1]). Here 
the use of a discrete Fourier transform D ; 91 = C[G] —91' = (S)iLi and 

the transfer of ideal decomposition problems to 91' is the most important tool to 
surmount difficulties. Calculations in 91' have the following advantages: 

(1) Decompositions (|4.1|) and solutions of (^4.2|) can be constructed very fast in 
91'. 


( 2 ) 

(3) 

(4) 


Every product formed during a run of (L) or (R) contains a factor which is 
a primitive idempotent e G 91'. Since every such e has only 1 non-vanishing 


block matrix E E 


i.e. 


e = (0,..., 0, E, 0,..., 0), the costs for every 


step of (L) or (R) reduce to the costs of calculations in a ring 
The algorithms (L) and (R) can be carried out completely within 91'. Only 
input and output data have to be mapped between 91 and 91' by means of 
D and D~^. Thus, ’’less fast” Fourier transforms can be useful, too. (See 
B. Fiedler 


lOi Sec.1.5.1].) 


In 91' there is a fast construction of bases of linear subspaces W, which we 
need to form linear equation systems 


To describe this construction, we denote by Ci^a ^ that matrix in which the 

Tth row is equal to a given a E K"" whereas all other rows are hlled with 0. 


Proposition 5.1. Let 1 = ■ A be a minimal left ideal o/ with known 

generating element 0 ^ A E and B = [5p]n,n 0 be a matrix from 

Determine a row a ^ 0 of A and a parametric form A of the solution of the linear 
eguation system 

n 

(5-1) E bijXj = 0 , i = 1,.. . ,n , Aj G IK (unknowns) . 
i=i 

Then B := {B ■ Ci^a \ 'I index for which A* is not a parameter in A} is a basis of 
the K-vector space B ■ i = B ■ • A. 


See B. Fiedler 


l0| . Sec.1.1.2] for the proof and other fast basis constructions. 
Spaces with a structure B ■ 1 are typical examples of spaces W (see Sec. i 0)- 
Further, we see that 1 has the basis B = {Ci^a | i = 1,..., u} if we use the identity 
matrix R = Id G for B. 

For our tensor investigations we need 91 = K[iSr]- M. Clausen und U. Baum 
j^, ^ developed a very fast Fourier transform for ]K[5j,], which bases on Young’s 
seminormal representation of Sr (see also H. Boerner |p and A. Kerber 1|T7| , Vol.I, 
p.75,76]). However, since the interpreter Mathematica does not allow the full speed 
and the optimal storage handling of this ingenious algorithm, we use Young’s nat¬ 
ural representation of Sr as discrete Fourier transform in our Mathematica package 
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PERMS 1^. (See H. Boerner pp. 102-108], B. Fiedler |Ty, Sec.1.5.2].) This 


implementation works good at least for Sr,r < 8. 

Multiplicities. Obviously, the efficiency of the algorithms (L) and (R) can be 
improved if we know before a run of (L) or (R) the multiplicities of equivalent 
minimal left/right ideals 1* or r* within decompositions 1 = 0™]^ 1* or r = 0™,]^ T 
searched for. If the algorithms have constructed such a direct sum of minimal 
left/right ideals of a fixed equivalence class that the number of summands equals 
the known multiplicity for this class, then the investigation of the remaining ideals 
of the class can be cancelled. This reduces the calculation time. 

In the case of fH = ]K[iSr] such multiplicies can be calculated by means of the 
irreducible characters of Sr, Frohenius reciprocity, the Littlewood-Richardson rule 
and plethysms. The determination of the irreducible characters of Sr is possible 
by the Murnaghan-Nakayama formula. We implemented all these tools in our 
Mathematica package PERMS (see B. Fiedler |^, Sec. II.3-II.6] for descriptions 
of implementations). For plethysms we use a very efficient method of F. Sanger 
P^, pp. 29-33]. (See B. Fiedler |^, Sec.II.6.3].) 


6. Characterizing Left Ideals of Tensor Products 

We continue to list linear subspaces W C ]K[iSr] describing tensor symmetries. In 
the case of tensor products two types of products can be considered: rd)(g)_ _ 
with possibly different and T 0 ... 0 T. 


Proposition 6.1. Q Let b C ]K[iSr.J (z = 1,... ,m) be left ideals and Pdl g C 
Tr^V be ri-times covariant tensors from the symmetry classes characterized by the 
U. Consider the product 

(6.1) T := 0 ... (g) G TW , r := n + ... + . 


For every i we define an embedding 

( 6 . 2 ) ii : Sri , {Ls){k) := 


Ai + s{k — Aj) if rj_i < k <ri 
k else 


where Aj := tq + ... + ri_i and ro := 0. Then the Tf, of the tensor ( 1 ^) fulfil 

(6.3) V6 G W : Tfc G 1 := K[5,] ■ £{Ii ■ ... • L} = K[5,] ■ (li ® ® L) 

where b := q(b) are the embeddings of the b into ]K[iSr] induced by the q. If 
dim V > r, then the above left ideal 1 is generated by all G ]K[iSr] which are 
formed from tensor products (| 6 . 1 |) of arbitrary tensors ^ 7 ^*, 


^^See B. Fiedler Sec.III.3.2] and B. Fiedler [Q. 
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Proposition 6.2. Q Let [q ^ be a left ideal and T G Ti* C T^V be a tensor 

of order m from the symmetry class T[*. Consider the product 

(6.4) f := T® ...®T e TmnV. 


Then all Tb, b G i/™", He in the left ideal 


(6.5) [ := K[S^n] ■ C{h ■ ... ■ in ■ I'} = • (li (8) ... ® ® i') 


where I* := ij(lo) are embeddings of [q into ]K[5mn] which are formed by means 
of mappings (|6.2|) with ri = ... = = m and r = mn. Further [' denotes the 

1-dimensional ideal I' := C,{Yfq^Q q} ofKlQ] where Q C Smn is the subgroup 


( 6 . 6 ) 


Q := 



k-m—l \ 
s{k)‘rrL—l/ 


l<k<n 


G iSn 


G 


rsj 


<Sn ■ 


If dim V > m ■ n, then the above left ideal I is generated by all Tb G K[iSmn] which 
are formed from tensor products (|6.4| ) of arbitrary tensors T ^ T[*. 


Let ug ■ G ^ GL(K[G]) denote the regular representation of a finite group G 
defined by Wg(/) := g ■ f, g ^ G, f ^ ]K[G]. If we use the above left ideals b, 
[q, t to dehne subrepresentations a* := a := /3 := a’ 5 ^|[, then the 

representation [3 is equivalent to a Littlewood-Richardson product or a plethysmQ, 
respectively (see B. Fiedler |^, Sec.III.3.2]): 


(6.7) 

( 6 . 8 ) 


[ according to ( |6.3| ) 
[ according to (|6.5|) 


/3 

/3 


Oi #.. 

a © \n\ 


# T 


These results correspond to statements of S.A. Fulling et ah H- (ED and (|6.8|) 
yield valuable information about multiplicities if one wishes to apply the algorithm 
(L) to 1. 


7. Subspaces Characterizing Tensors with Index Contractions 

First we give a universal linear subspace which contains the group ring elements 
^ tensor T with I index contractions for every value of dim V. 

Theorem 7.1. Let V,B,r,l, g,bo have the meaning given in Def. and Prop. 
1^. Consider the partition Aq := ( 2 *, 1 '’“^^) h r and the lexicographically smallest 
standard tableau t of Xq. Form the group ^ G := Ht ■ Q where Tit is the group 
of all horizontal permutations of t and Q (Z Vt is the subgroup of all such vertical 

^^See B. Fiedler Sec.III.3.2] and B. Fiedler [p^. 

^^See the references |^, |4[ |lj for the Littlewood-Richardson rule and 

plethysms. 

is a semidirect product Ht Q and isomorphic to the wreath product S 2 IS 1 . 
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permutations of t which only permute full rows of t with length 2. Then every 
tensor T & T^* %¥ ({ = K[5r] • e, e idempotent) fulfils 


(7.1) 


^ e 1 g ■ IK[iSr] • e , 1 g ■— 9 ■ 


be'Bt.n 


g&G 


Furthermore, if dim V >r — I, then there is such a bo E B’' that 


(7.2) 


1g ■ IK[5r] • e — £]k| IbTb T G 




The proof can be found in Sec. III.3.4 of our Habilitationsschrift If dim V < 
r — I, then the Xlbe'Bi, 7^ span only a linear subspace of 1 g ■ IK[iSr] • e in gen¬ 

eral. To describe this subspace, we define: 


Definition 7.2. If A = (Ai,...,Afc) h r is a partition with length |A| = k and 
(ui,..., Ufc) G is a /c-tuple of vectors, then we denote by (A; Ui,..., Vk) or short 
(A; Vi) that r-tuple from which has the structure 

(7.3) (A;ni,... ,nfc) := ( ui,. , Ui , Ua,. ,Us ,..., , Vk) G V’". 

Ai A 2 Afc 

For every b = (vi,...,Vr) G V^, there exists a unique partition A h r and a 
permutation q E Sr such that b can be written as 6 = q{X;wi,... ,w\x\) where 
Wi,..., tC|A| are the pairwise different, suitably renumbered vectors from b. We call 
{X;wi,... ,tC|A|) a grouping of b and A the grouping partition of b, which we also 
denote by A = 6 '^. 

Definition 7.3. Let B be an orthonormal basis with respect to a fundamental 
tensor g E T 2 V. We call (uj^,... ,nj^,) G B^ smaller than (n^^,... ,nj^,) G B^ if 
the first non-vanishing difference jk — h fulfils jk — ik > 0. If (A; wi ,..., tC|A|) and 
(A; ..., are two groupings of a fixed r-tuple b E B^ of basis vectors, then we 

call (A; wi ,..., tC|A|) smaller than (A; w [,..., if the |A|-tuple (uj^,..., := 

{wi,... ,W|A|) is smaller than the |A|-tuple (n^^,... := {w[,... ,u'jA|)- 

For every r-tuple b E B^ there exists a permutation pb E Sr such that b has a rep¬ 
resentation b = Pb{X; Wi, ..., w\\i) where {X;wi,..., tCjAi) is the smallest grouping of 
b and X = b^. We denote by p a single-valued mapping p : B^ ^ Sr,b ^ p{b) '. = Pb 
which assigns exactly one of such permutations pb to b. 

Let bo E B^~^'' be an (r — 2/)-tuple of vectors from the basis B. We denote by 
AbQ the set Ab^ := {A h r | 36 G fSbo ■ ^ Furthermore, we assign to every 

partition A G the lexicographically smallest standard tableau tx of A and the 

^^Corollary: If r — 2Z = 0, then the decomposition of Iq ■ IK[5r.] into minimal right ideals is 
characterized by a plethysm [2] © [1] ^ Thus the number I of linearly independent 

invariants of T is bounded by the sum M of the multiplicities of minimal left ideals b belonging 
to partitions 2g, g E I, in IK[5r.] • e = 0b. If dimF > I, then / = M. (See B. Fiedler 0 
Sec.III.4.2]. Compare S.A. Fulling et al. Q.) 
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set A4bo,\ ■= {p(^) E S'" \ b E ^bo with b^ = A} of such r-tuples which are the 
smallest groupings of the b E ^bo with grouping partition A. 


Theorem 7.4. Let g,bo,T[* and e have the meaning given in Theorem 


7J\ and p : ^ Sr be a mapping of the type described in Definition Then we 


have 


7b ^6 T G 7^*1 — ^ ^ «(A;u>i> ■ Iwt;, ■ • e 


b&<Bb^ 


where 


a{x,u,i) ■= ^bp{b) K 

bG®6o 

p{b)-^b= (X]Wi) 


The proof is given in Sec. III.3.4 of our Habilitationsschrift [|T^. Furthermore, 


we 


clear in Sec. III.3.4 to what extent the subspaces from Theorem ^]Jare independent 
from the choice of bo. 


8. Concluding Remarks 


rm we determined linear subspaces 


12 , too. 


10] we applied Theorem 7.1 to invari 


In Sec. III.3.3 of our Habilitationsschrift 
W C ]K[i5r] for tensors T G TrV for which the vector space V has a dimension 
dim V < r. Some details about this case can be found in the paper 
In Sec. III.4.2 of the Habilitationsschrift 
ants (i.e. r = 21 and 6o = 0). Footnote |T^ gives one of the results. 

We have tested our methods in computer calculations. Among other things, we 
treated the term combination problem for quadratic monomials in the coordinates 
of the Riemannian curvature tensor (B. Fiedler ||I0| , Sec.HI.5.1]). The correct¬ 
ness of this calculation was controlled by means of the results of S.A. Fulling 
et ah 1^. Furthermore, we used Theorem o to verify the standard identitfi^ 
J2p^S2n a(7')^p(i) ■ • • • • Ap( 2 n) = 0 for (n x n)-matrices Aj in the case n = 2 (B. 
Fiedler [^, Sec.HI.5.2]). For this calculation the full algorithm (R) was required. 

We implemented our methods in a Mathematica package PERMS. This package 
comprises tools for permutation groups and group rings, partitions and tableaux, 
special idempotents, characters, the Littlewood-Richardson rule, plethysms, dis¬ 
crete Fourier transforms and our algorithms (L) and (R). Furthermore, we used 
the packages SYMMETRICA ]]T^ and GAP 
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